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ABSTRACT 

Carbon monoxide (CO) is the most commonly used tracer of molecular gas in the inner regions of 
protoplanetary disks. CO can be used to constrain the excitation and structure of the circumstellar 
environment. Absorption line spectroscopy provides an accurate assessment of a single line-of-sight 
through the protoplanetary disk system, giving more straightforward estimates of column densities 
and temperatures than CO and molecular hydrogen (H2) emission line studies. We analyze new 
observations of ultraviolet CO absorption from the Hubble Space Telescope along the sightlines to 
six classical T Tauri stars. Gas velocities consistent with the stellar velocities, combined with the 
moderate-to-high disk inclinations, argue against the absorbing CO gas originating in a fast-moving 
disk wind. We conclude that the far-ultraviolet observations provide a direct measure of the disk 
atmosphere or possibly a slow disk wind. The CO absorption lines are reproduced by model spectra 
with column densities in the range N{^'^CO) - 10^^ - 10^® cm^^ and N{^^CO) - lO^^ - 10^^ cm'^, 
rotational temperatures Trot{CO) ^ 300 - 700 K, and Doppler 6-values, 5 ~ 0.5 - 1.5 km . We 
use these results to constrain the line-of-sight density of the warm molecular gas [nco 70 — 4000 
cm~^) and put these observations in context with protoplanetary disk models. 

Subject headings: protoplanetary disks stars: individual (AA Tau, HN Tau, DE Tau, RECX-15, 

RW Aur, SU Aur) 



1. INTRODUCTION 

Characterizing the gas properties of protoplanetary 
disks at planet-forming radii (a < 10 AU) is criti- 
cal to understanding the way planets form and evolve. 
The lifetimes of protoplanetary disks are comparable 
to the time during which giant planet cor es form and 
accre te gaseous envelopes ( 10^ — 10^ vr: iFedele et all 
IMoHHer^andez et al.l 120071 : iHubickvi et al.ll2005D . The 
surface density of the disk in which the planet forms 
plac es limits on the amount of migration tha t is possi- 
ble ()Trilling et al.|[200llAmiitage et al.|[200l . The life- 
time of the protoplanetary disk determines the amount 
of time that giant planets are able to accrete, thereby in- 
fluencing their final masses (.Ida fc Lin .2004 ). Observa- 
tions of carbon monoxide (CO) emission and absorption 
lines from ultraviolet (UV) to infrared (IR) wavelengths 
have been shown to be useful probe s of molecul ar gas 
in the inner disk (Naii ta et all 120031: ISalvk et al.l 120111: 
IKchindhclm ct al. 2012). 

Emission and absorption from the CO Fourth Posi- 
tive band system {A ^11 — AT ^S+) has been widely used 
to study the inters tellar medium (iFederman et al.lll980l: 
Burgh et all |200/), debris disks ("Vidal-Madia r et al ' 



1994 ). comets (Fcldman & Brunc 1976; McPhat e et al 
19971 ). planets (lDurrancelll981l : IFeldman ct al. 200Ci), and 



the atmospheres of cool stars (jCarpcntcr et al.. il994[ ). 
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The first detections of CO far-UV emission and absorp- 
tion spectral features from the i nner regions of proto - 
planetary disks were presented in iFrance et al.l (|2011aD . 
The detected CO absorption lines have rotation temper- 
atures Trot{CO) PS 500 ± 200 K, pointing to an origin 
in the warm inner disk gas. Molecular hydrogen (H2) 
and CO absorp tion were me a sured simultaneously for 
the first time by IF rance et al.l (j2012aD in the disk of AA 
Tauri. The H2 w as seen in absorpti on against the Lya 
emission line, as in lYang et al.l (j2011l ). and the CO A — A 
absorpti on bands were obser ved against the far-UV con- 
tinuum (IFrance et al.ll2011bD. The value of CO/H2 - 0.4 
found bv lFrance et al.l (|2012aD is approximately three or- 
ders of magnitude larger than the canonical interstellar 
value of 10""*. Better constraints on the CO/H2 ratio 
are important for determining the total amount of gas in 
protoplanetary disks. 

Many classical T Tauri star (CTTS) disks have been 
studied with near-IR fundamental (4.7 ^va) and over- 
tone (2.3 ftm) bands of CO, in both emission (e.g. 
ICarr et al.l 11991 [Najita ct al. 2003; Salvk ct al. 201^ 
and absorption (e.g. IRettig et al.i l2006l : IHorne et all 
120121) . Near-IR CO emission has been used as a tracer 
of the dominant molecular gas component, H2, because 
the homonuclear nature of H2 makes rovibrational tran- 
sitions dipole-forbidden and thus difficu lt to observe 
at near- and mid- I R wavelength s ( Pascu cci et al.l 120061 : 



Bitner et al.l 120081; ICarm ona et al. 2008i: althoug h see 
2003; Ramsavl fo^at fc Greavei [2007: and 
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The CO lines at 4.7 /im appear to 
originate from the inner disk region with T « 500 - 
1500 K and Keplerian velocities consistent with the dust 
sublimation radius. However, some source s have line 
profiles with excess low velocity emission (|Bast et al.l 
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[20Tl[) . which iPontoppidan et all (|20Tll) propose arises 
from a slow molecular disk wind. In general, the IR 
CO temperatures fr om LTE slab mod els (275 - 1675 K, 
average ~ 1100 K; ISalvk et all [20Tll) are higher than 
those measured from UV CO fluorescence (460 ± 250 
K: ISchindhelm et al.l [2011. iRettig et all (|2006l ) and 
iHorne et al.l ( 20121) studied IR CO absorption from the 
surfaces of circumstellar disks. IRettig et al.l (|2006f ) find 
CO excitation temperatures of ~ 100 K, likely placing 
the absorbing gas farther out in the disk at r > 10 AU. 

The discrepancies between the IR and UV CO work 
suggests that the observations are probing multiple 
molecular gas populations, and multiple populations may 
be responsible for the anomalously high observed CO/II2 
ratio in AA Tau. With deep UV absorption spectroscopy, 
we have a new observational tool for studies of the inner 
disk gas. Absorption lines give the most direct measure 
of the column density and temperature of the disk gas 
along the line-of-sight as the line fitting is largely inde- 
pendent of the geometry or the photo-exciting source. In 
order to provide observational constraints on inner disk 
gas, we present new analyses of far-UV spectroscopic ob- 
servations of six (0.6 - 6 Myr) CTTSs. Our observations 
probe the warm disk atmosphere, which is important for 
constraining the three-dimensional structure of disks in 
which planets form. We describe the targets and the ob- 
servations in §2. The analysis performed and the model 
fit parameters are presented in §3. The disk geometry 
and correlations in the data are discussed in §4. Finally, 
§5 contains a brief summary of our results. 

2. TARGETS AND OBSERVATIONS 

We analyze far-ultraviolet spectra of six targets that 
are a subse t of th e 34 T Tauri stars presented in 
[France et all (|2012b[ ): AA Tau, HN Tau A, DE Tau, 
RECX-15 (ET Cha), RW Aur A, and SU Aur. HN Tau 
A and RW Aur A (hereafter HN Tau and RW Aur) are 
binary stars, but with separations large e nough, 1.4" for 
RW Aur AB and 3.1" for HN Tau AB (|White fc Ghed 
1200 It ), such that only the primary is within the aper- 
ture and dominates the emission analyzed here. AA Tau 
and HN Tau have been studied previously by France et 
al. (2011a, 2012a), but we refit their spectra with the 
same procedure as the other targets for consistency. Five 
targets are in the Ta urus- Auriga star-forming region at 
a distance of 140 pc (jBertout et al.l [19991 iLoinard et al.l 
[M07). RECX-15 is in the 77 Cham aeleontis cluster at 
a distance of 97 pc (jMamaiek et al.1 '1999'). Parameters 
including the age, luminosity, and spectral type of these 
targets are listed in Table 1. 

All targets were observed using the Hubble Space 
Telescope-Cosmic Origins Spectrograph {HST-COS; 
[Green et al.]|2012[ ) under program ID 11616. The COS 
FUV M-mode wavelength solution is accurate to Av ^ 
15 km s~^ and depends on the object centering. Target 
acquisition was through the MIRRORB near-UV imag- 
ing mode for AA Tau and DE Tau. The rest of the tar- 
gets were acquired using the MIRRORA mode. Far-UV 
spectra were obtained using three central wavelengths 
for G160M, two central wavelengths for G130M, and 
multiple focal-plane positions to cover the 1133 < A < 
1795 A bandpass while minimizing fixed pattern noise. 



' http: / / www.stsci.edu/hst / cos / ducuments /handbooks / current / cos_cover 



Table 2 lists the dates and exposure times for each ob- 
ject used in this paper. The data were processed through 
CALCOS, the COS calibration pipeline, and aligned and 
CO- add ed with the procedure described in lDanforth et all 

The six targets were the only objects in the T Tauri 
star sample to show unambiguous CO absorption. The 
rest of the larger sample either had low continuum signal- 
to-noise (S/N) making modeling difficult, or were com- 
plete non-detections of CO absorption (most likely from 
low-inclination disks not intercepting the line-of-sight). 
Our optical depth detection limit was r = 5, below which 
the absorption depths were comparable to the noise in 
the data. The targets that we do not fit have incli- 
natio ns ranging from 4° (TW Hya, I Pontoppidan et al.l 
[2OOI to 85° (DF Tau, |JolmsJirulL& Valcnti 2001) with 
an average of 43°. Thus, the majority of the targets 
that do not show strong CO absorption have lower in- 
clinations. DF Tau is an unusual target discussed in 
§4.3. The six targets that we analyze all have medium- 
to-high inclination (35° - 77°, with an average of 61°) 
protoplanetary disks. The average does not include HN 
Tau, whose inclination is not well known. The inclina- 
tions have been determined in a number of different ways, 
with fairly large uncertainties arising from the techniques 
used. The inclination of AA Tau was determined by pe- 
riodic eclipses of the star by the magnetically warped 
accretion disk combined with a measure d line-of-sight 
projected rotation velocit y in the optical ([Donati et al.l 
[20Ta IBouvier et al[[2003[) . For DE Tau, it was calcu- 
lated from stellar radii and rotation periods spectroscop- 
ically determined in the I R combined with the lite rature 
value of V sin i in the red pohns-Kruh fc Valenti| [2001'). 
Optical spectroscopy of a very strong H^ emission line 
profile was m odeled to find the inclination of the disk 
of RECX-15 ([Lawson et al.l 120041 ). The interferometric 
data of RW Aur was fit with an inclined uniform disk 
model an d was sufficient to constrain the inclination of 
the disk (jE isner et al.ll2007[ ). This inclination IS mcon- 
sistent with the RW Aur inclination 45°) found by 
iLopez-Martfn et al.l ([2003| ) using the ratio of proper mo- 
tion to radial velocity toward emission knots in the RW 
Aur jet. The visibility as a function of hour angle in the 
K-band of SU Aur was fit with a Gaussian brightness 
profile that was inclin ed on the sky, givin g the best fit 
inclination for the disk (Ake son et al.l[2002[ ). The inclina- 
tion of HN Tau is not well known, and only a lower limit 
is given to indicate that an inclination higher than 40° 
is needed for the line-of-sight to intercept the di sk and 
to explain observed outflows ([France et al.ll2011a[) . This 
lower limit is increased in §4.3 by comparing to published 
models. 

The radial velocities of the majority of the stars (listed 
in Table 3) were obtained bv lNguven et al.l (|2012l ). The 
radial velocit y for RW Aur A was not measured by 
[Nguven et al.l ([2012D . but is assumed to be the same as 
that for the secondary, RW Aur B. The radial velocity of 
RECX-15 is determined using photospheric Li I absorp- 
tion lines at 6707.76 A and 6707.91 A ([Woitke et al.l 
1201 ID with resolving power of - 30, 000 (10 km s'^ at 
6300 A). 

The CO is observed in absorption against the far- 
UV, continuum of these actively accreting stars. The 
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TABLE 1 

Target Parameters 



Object 


Spectral 




Inclination 


L, 


M. 


M 


logi()(Age) 


Ref.^' 




Type 




(degrees) 


(L©) 


(M©) 


(10-8 M0 yr-i) 


(yrs) 




AA Tau 


K7 


0.5 


70 


0.71 


0.80 


0.33 


6.38 ± 0.20 


2,4,7,12,16 


DE Tau 


MO 


0.6 


35 


0.87 


0.59 


2.64 


5.82 ± 0.20 


7,10,12 


HN Tau 


K5 


0.5 


>40 


0.19 


0.85 


0.13 


6.27 ± 0.27 


6,7,12 


RECX-15 


M2 


0.0 


60 


0.08 


0.40 


0.10 


6.78 ± 0.08 


13,14,15 


RW Aur 


K4 


1.6 


77 


2.3 


1.40 


3.16 


5.85 ± 0.53 


5,9,11,12,17 


SU Aur 


Gl 


0.9 


62 


9.6 


2.30 


0.45 


6.39 ± 0.21 


1,3,8,11,12 



^ ri'l lAkeson et al.l (I2002D : (2) [Andrews fc Williamj (I2007D: (31 iBertout et al.l (I1988D: (41 iDonati et all (l20ldF 
(5) lEisnereFall pOOl : ( 61 [France et al.l (1201 lal); (7) IGullbring et"a r (lm%: (81 [ GuUbr ing et al.l (l2000ll: (91 
IHartiean et all (1199511: (101 'Johns-Krull fc Valentil (|200l]l: (l ll lJoWKrull et all mm~(l2) IKraus fc Hillenbrandl 
([2f)09h: (131 ILawson et al. (2004); (14) Lu hman fc Steeghsl (pOOl : (15) IR.amsav Howat fc Greave^ (pOOl : (16) 
IRicci et alJ (|20Toll : (17) IWhite fc Ghea (1200111 . 



TABLE 2 

HST-COS G130M/G160M Observations of Targets 



Object 


R. A. (J2000) 


Dec. (J2000) 


Date 


G130M Exposure 
(s) 


G160M Exposure 
(s) 


AA Tau 


04 34 55.42 


+24 28 52.8 


2011 Jan 06, 07 


5588 


4192 


DE Tau 


04 21 55.69 


+27 55 06.1 


2010 Aug 20 


2088 


1851 


HN Tau 


04 33 39.37 


+17 51 52.1 


2010 Feb 10 


5724 


4528 


RECX-15 


08 43 18.43 


-79 05 17.7 


2011 Feb 05 


3890 


4501 


RW Aur 


05 07 49.51 


+30 24 04.8 


2011 Mar 25 


1764 


1617 


SU Aur 


04 55 59.39 


+30 34 01.2 


2011 Mar 25 


1788 


1759 



far-UV continuum arises from stellar photospheric and 
chromos pheric emission as • well as an accretion con- 
tinuum (iCalvet fc GullbTi^ [19981 iFrance et al]|2011bl: 
llnglebv et al.ll2011l : see Table 1 for the accretion rates). 
The far-UV continuum may be enhanced in some 
objects by a m olecular dissociation quasi-continuum 
(jBergin et ani2004l: iHerczeg et al.ll2004l ). but the details 
of this process have not been thoroughly ch aracterized 
observationally (see, e.g.. lFrance et al.ll2011b( ). The con- 
tinuum emission is absorbed by circumstellar atoms and 
molecules. Depending on the inclination angle, the light 
may pass through multiple layers of the disk and thus 
probe multiple temperature and density regimes. At the 
high temperatures and column densities of the observed 
CO (see §3.3), the contribution from interstellar CO is 
negligible for transitions from J > 2. For temperatures 
typical of the interstellar medium (ISM), states with J 
> 3 are not significantly populated. For our extinc- 
tion values (Ay ^ 1-5), we would expect an interstellar 
CO column densi ty of 7V(C0) < 10" (see Figure 3 of 
iBurgh et al.ll2007l ). which is orders of magnitude smaller 
than the circumstellar column densities that we observe. 
Therefore, our CO measurements do not probe the cold 
diffuse interstellar material with <Trot{CO)>isM ~ 4 K 
(|Burgh et al.ll2007D . but arise from material close to the 
star, likely the warm molecular disk. 

3. ANALYSIS 
3.1. Data Analysis 

We identify nine CO bands, {v - 0), that span the COS 
far-UV bandpass from the (0 - 0) band at A ~ 1544 A to 
the (8 - 0) band at A ~ 1322 A (Figure 1). We do not 



fit the (5 - 0) band due to contaminating emission from 
H2 and Si IV A1394 A. The (0 - 0) band is contami- 
nated with C IV A1548 A emission and is also excluded 
from the overall fit. Many photo-excited H2 emission 
lines populate the spectra, complicating the model spec- 
trum fits as well. All CTTS systems with identifiable 
CO absorption contain broad bands created by the over- 
lap of many closely spaced rotational lines. This sug- 
gests a relatively high rotational temperature, Trot{CO), 
because as higher-J lines are populated, the observed 
spectral bands become broader. Broad CO absorption 
bands were first ide ntified in the spectrum of HN Tau 
(iFrance et al.l 1201 lal) and later in the spectrum of AA 
Tau ( France et al.ll2012al ). We employ a similar tech- 
nique to the previous work, using spectral synthesis mod- 
cling to measure Trot{CO) and N(CO) in all six sources. 

The spectrum around each band was normalized by a 
linear fit to the nearby continuum, avoiding regions with 
II2 or other emission lines. The least contaminated ab- 
sorption bands are the (1 - 0), (2 - 0), (3 - 0), and (4 - 0) 
bands, which also have the strongest oscillator strengths 
in the Fourth Positive system. However, ^^CO and ^-^00 
are not cleanly resolved in the low-w bands. To isolate 
these isotopologues, higher-w bands with larger separa- 
tions between the ^^CO and ^'^CO bandheads (f = 7, 
8) are used to better constrain the column densities of 
both species (see §3.3). There is nearly a 3 A separa- 
tion between the isotopes in the (7 - 0) band, compared 
to the '^1 A separation in the (1 - 0) band. However, 
these higher bands can sometimes blend into the contin- 
uum because of their shallower absorptions. HN Tau in 
particular does not have detected w = 7, 8 absorption 
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Fig. 1. — The nine far-UV CO rovibrational absorption bandheads identified by orange vertical fines. Going from left to right: {v - v 
= 8 - 0) 1322.1 A, (7 - 0) 1344.2 A, (6 - 0) 1367.6 A, (5 - 0) 1392.5 A, (4 - 0) 1419.0 A, (3 - 0) 1447.4 A, (2 - 0) 1477.6 A, (1 - 0) 1509.7 A, 
(0 - 0) 1544.4 A. The spectra have been Fourier smoothed to suppress frequencies higher than the 7-pixel spectral resolution. 
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Fig. 2. — how-v (above) and high-ii (below) model fits for 
DE Tau showing the larger separation between ^■^CO and ^^CO 
for larger v -values. The separation helps to better constrain the 
^■^CO/ CO ratio in the model. The data is in black and the best- 
fit model is in orange. The first absorption dip in each model is the 
^^CO bandhead, and the ^^CO bandhead is marked with a blue 
arrow. 

(see §3.3). An example of these higher band fits showing 
separations in the DE Tau (1-0) and (7-0) fines is 
sfiown in Figure 2. 

Tfie normafized data vary in both the depth of their 
absorptions and the number of excited J- fines present. 
DE Tau and SU Aur have the largest depths, but the 
data are also noisier, making it harder to determine a 
continuum level. The depths will vary, even for similar 
column densities, because of the error associated with 
the continuum fit and population shifts in the lower- J 
levels due to increases and decreases in temperature (see 
§3.2). Because the absorptions in these targets are very 
saturated (see §3.3), changing the column density does 
not affect the absorption depth highly at low- J values. 
Here, the normalization will have a bigger effect. At 
high- J values (longer wavelengths within the band) it is 
difficult to pinpoint tfie weaker absorptions amongst tfie 
emission lines and tfie noise in the data. This makes 
the data harder to interpret, though RW Aur does have 
a significant lack of high- J lines compared to tfie otfier 
targets. HN Tau, however, has fairly prominent high- 
J lines. All of these factors contribute to the derived 
molecular parameters and errors described in §3.2 and 
§3.3. 

3.2. Model Description and Fitting Procedure 

The CO absorption profiles were modeled for ^^CO 
using the oscillator strengths and rovibrational line 
wavelengtfis from Eidelsberg (private communication). 
The A energy levels of iHaridass fc Hubed ([1991 
were used to calculate the ground state energies. 
The oscillator s treng ths for ^'^CO are taken from 
lEidelsberg et al.l (|1999l ). and the wavelengths and oscil- 
lator strengths_Joi_Jhfi_2erto states were taken 
from lEidelsberg fc Rostai (|200l . A common rotational 
temperature was assumed for both ^^CO and ^"^CO. 



All six of the analyzed bands (u = 1, 2, 3, 4, 7, 8) were 
compared to model spectra for a grid of different values of 
the Doppler 6- value, rotational temperature, Trot{CO), 
and logarithmic column density of both -"^^CO and ^'^CO, 
\ogw{N{^^GO)) and logio(iV("CO)). The model also in- 
cludes a velocity shift parameter (explored in more de- 
tail in §4.1), which we originally estimated by eye from 
the (2 - 0) band of each target by varying the velocity 
shift until the bandheads of the model and data coin- 
cided. However, after a more rigorous treatment of the 
velocity shifts, using all the bands with uncontaminated 
bandheads to constrain the location of the CO gas in the 
disk, more reliable values are found in §4.1. The ranges 
of our grid search were 0.5 - 2.0 km s~^ in steps of 0.1 km 
s~^ for the Doppler 6- value, 300 - 1000 K with steps of 
50 K for the rotational temperature, 16.0 - 18.0 in steps 
of 0.1 for logio(A^(i2CO)), and 14.0 - 17.0 in steps of 0.1 
for logio(Af(^^CO)). The maximum 6-value for our grid 
search comes from our assumption of turbulent velocities 
< 1 km s~^ and CO rotational temperatures < 5 x 10^ 
K. The maximum login(i V(^^CO)) value for t he search 
comes from the analysis of lFrance et al.l (j2012aD , who ar- 
gue that any larger column of ^^CO would lead to a col- 
umn of hydrogen that would extinguish the stellar flux 
around the Lyman alpha line center, wfiicfi is inconsis- 
tent witfi our observations. This limit is confirmed after 
the fit as tfie majority of tfie targets (except those with 
the highest column, RECX-15 and SU Aur) are incon- 
sistent with a column of 10^^ cm~^ within their errors. 
Any ^^CO column density lower than our lower limit of 
10^^ cm^^ (wfiicfi corresponds with our optical depth de- 
tection limit of r = 5) would be very difficult to detect 
witfi tfie S/N, setting the lower bound on our column 
density search. An example of our CO model is shown 
in Figure 3. The higher rotational states are excited at 
higher temperatures while the lower rotational states are 
depopulated. Thus, the observed absorption band be- 
comes wider and shallower with increasing temperature 
if column density is kept constant. The three simulated 
spectra at the top of Figure 3 are at the native model 
resolution (Av ~ 2 km s~^) and the bottom spectrum 
is at the resolution of the HST-COS instrument (Av ^ 
17 km s~^). The CO absorption lines are narrower than 
the instrumental resolution. Tfie unresolved line cores do 
not go to zero at the observational resolution even when 
the CO optical depths are large. 

Because the statistical errors in the data were some- 
times anomalously smalQ, such that a reduced fitting 
routine with non-uniform errors gave erroneous results, 
we employed a simple least squares fitting routine that 
minimizes the square of the difference between the model 
and the data. This procedure weighted all data points 
equally, giving a more reliable fit. Because of the pres- 
ence of emission lines (mostly H2), we truncated the data 
above a normalized flux level of 1.1 so that tfie emis- 
sion features did not contaminate the absorption model 
fitting. This does, however, cause the model to fit a 
truncated continuum value instead of tfie absorption at 
times when the emission lines wipe out the absorption. 
The model is pulled toward a solution with lower column 

^ The CALCOS pipeline is known to mishandle statistical errors 
in the low S/N regime iFroning et al.|[2011 ) 
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Fig. 3. — Three different CO models for tlie (2 - 0) band 
of DE Tau with identical Doppler 6-value = 1.4 km s~^, 
logio(A^(^^CO)) = 17.2 , logio(Af("CO)) = 16.0, and velocity 
shift = 1.4 km . (Top 3): Model output with native model 
resolution (Av ~ 2 km s~^). (Bottom): Model output after being 
convolved with the HST-COS linespread function (Av ~ 17 km 
s~^) to simulate data detected with the instrument. 



density to fit the higher flux values at these points, in- 
troducing a bias which contributes to the large errors on 
the column density in Table 3. 

The model uses the five input parameters to create a 
theoretical spectrum. The Doppler 6-value, along with 
the column density, determines the opacity of the model, 
and the temperature determines the band shape as dis- 
cussed above, with higher temperatures leading to wider, 
shallower absorption models. The column density inputs 
determine the relative depth of the ^^CO and ^^CO ab- 
sorption bands. The ^^CO/^^CO ratio is primarily de- 
termined by the differences in the absorption depths at 
the relevant bandheads. At the high column densities of 
these disks (see §3.3), most of the targets are on the flat 
portion of the curve of growth (COG) for the lower- J 

^^CO lines of the lower-?; bands. This means that in- 
creasing the column density of ^^CO does not lead to 
a significant increase in the absorbed flux. The column 
density fits rely more on the "'^'^CO absorption and the 
higher-w bands and higher-J lines of ^^CO, which are 
less saturated. The velocity shift parameter moves the 
absorption band of the model to shorter or longer wave- 
lengths corresponding to the gas moving toward or away 
from the observer, respectively. The theoretical spec- 
trum is then convolved with the COS linespread function 
and compared to the observations to find the best-fit pa- 
rameter values. 

3.3. Model Fits and Errors 

The fits to the normalized absorption data are shown 
in Figure 4 and the parameter values are provided in 
Table 3. DE Tau and SU Aur, which have the deepest 
absorptions, have higher column densities in both ^^CO 
and ^^CO than most of the other targets. The shal- 
lower low- J lines combined with the prominent high- J 
lines in HN Tau require a higher temperature and lower 
column density to fit the data. Most of the targets 
share a common best-fit temperature (^ 350 K), with 
HN Tau being an outlier because of the more prominent 



1.5 
1.0 
0.5 

X 

I 0.0 

0) 
N 

S 

o 1.5 

Z 

1.0 

0.5 
0.0 




CO (1-0) 



i 

,1 0J 



CO (2-0) 



1509 1512 1515 1518 



1477 1480 1483 1486 



COi}-0) 



f \ i 



J' I' 



CO (4-0) 



1448 1451 1454 1457 1418 1421 1424 1427 
Wavelength (A) 

Fig. 4a. — Low-D model fits (orange) for AA Tau. The 
best-fit values (see Table 3) for the Doppler b-value, Trot{CO), 
logio(Ar(i2cO)) , logio(Ar("CO)), and velocity shift are used. The 
data (black) are continuum normalized and the ^^CO bandheads 
are marked to clearly show the two different CO species. The emis- 
sion contaminating the absorption is not included in the model and 
any flux above a normalized value of 1.1 is truncated during the 
fit. 
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Fig. 4b. — Same as Figure 4A for DE Tau. 



high- J lines. The column densities are much more var- 
ied, with large uncertainties because of the high- J lines 
blending in with the low S /N continuum and contamina- 
tion by emission lines. Unlike the temperature, there is 
no common best-fit column density for all of the targets 
within the uncertainties. An intersystem band not taken 
into account in the code (a' - X (14 - 0)) at 
A = 1419.50 A (jEidelsberg fc Rostasl[200l hinders the 
fit of the (4 - 0) band of HN Tau, the highest temper- 
ature target. This absorption is between the ^^CO and 
^■^CO bandheads, causing the data to be deeper than the 
model in this part of the spectrum. A few targets also 
show this difference in the (2 - 0) band as well, which 
could be another intersystem band. 

Our fit values for all six targets lie in the ranges 
of 6 = 0.5 - 1.4 km s-\ Trot{CO) = 300 ~ 700 K, 
logio(iV(i2cO)) = 16.0 - 17.8, and logio(Af("CO)) = 
14.7 - 16.7. The errors were estimated by taking the 
standard deviation of the distribution with appropri- 
ate degrees of freedom. Increasing our minimum value 
by the standard deviation determined the best-fit param- 
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TABLE 3 
CO Fit Parameters 



Object 



(10- 



t(1420) 



erg cm 



logio(Ar{i2cO)) logio(Af("CO)) 



TrotiCO) b vco 

(K) (km s-l) (km s-l) 



v» V* Reference'' 

(km s"-"-) 



AA Tau 


6.1 


16.9l°f 


15.5+1'° 


450ti0 


1 4+0'7 
-"-•^-0.6 


+24.7 ± 18.6 


16.98 + 0.04 


2 


DE Tau 


3.7 


17 9+0.6 


16.0+i-O 


450i^r 


1 4+0.6 


+ 1.4 ± 15.9 


15.402 + 0.018 


2 


HN Tau 


7.9 


16.01°:^ 


15.2+?-^ 


7nn+i2oo 

'UU_45o 


1 4+1-3 


+28.0 + 21.1 


4.6 ± 0.6 


2 


RECX-15 


4.5 


i7.8ig:^ 


16.71^0 


oc-n+1850 
■JOU_2oo 


U.O_(j 2 


+18.2 ± 11.3 


22 


1 


RW Aur 


22.3 


i6.9+;« 


14.71^7 




6+''-'' 

U.D_Q 2 


+20 ± 20 


15.00 + 0.03'' 


2 


SU Aur 


5.6 


17 4+0.7 


15.511:^ 




1-^-0.5 


+16.3 + 27.3 


14.26 + 0.05 


2 



fP IWoitke etaP |[20TlD : C2'l lNguven et al.l S20W \. 

^ Value for RW Aur B 
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Fig. 4c. — Same as Figure 4A for HN Tau. 
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Fig. 4e. — Same as Figure 4A for RW Aur. 
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Fig. 4d. — Same as Figure 4A for RECX-15. 
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Fig. 4f. — Same as Figure 4A for SU Aur. 



eter range. We define errors as the wi(ith of tfie best-fit 
parameter space, as defined by tfie standard deviation of 
the distribution. This error procedure is illustrated 
in Figure 5. The errors on logio(iV(^'^CO)) were more 
difficult to define because of the lower column density of 
13C0. For logio(iV(^^CO)) < 14, the value does not 
increase substantially as the column density is decreased 
because we are no longer detecting the absorption. We 
take a lower limit on the column density of ^^CO of 10^^ 
cm~^ for this reason. Within the errors, HN Tau, RW 
Aur, and SU Aur are consistent with a non-detection 
of i^CO (column density of "CO < lO^'* cm'^). The 



temperatures in general tend to have large upper error 
bars. Higher temperatures populate the higher J-states, 
which produce weaker absorption lines due to smaller 
column densities. Due to the low S/N in the data, these 
weak lines blend into the continuum making the high 
temperature limit difficult to constrain. At temperatures 
> 1000 K, the {v - 1) rovibrational bands of CO should 
be detectable (t > 5) in the data at our best-fit column 
densities. At our best-fit temperature of ~ 350 K, the 
{v' - 1) band would be detected for logio(7V(i2CO)) > 
18.0. However, we do not detect any of these bands. 
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Fig. 5. — Illustration of the uncertainty determination. (Above): 
A plot of the distribution for the degrees of freedom of RW 
Aur. The standard deviation of the distribution is shown in red; 
the mean is defined by the dashed line. (Below): A plot of the 
unweighted value as a function of logio(A'^(^^CO)). We increased 
the minimum x^ by the standard deviation of the distribution to 
define our error region. The dashed line shows the best-fit column 
density for RW Aur (10^^'^ cm~^), the red line shows the standard 
deviation of the above x^ distribution, and the orange line defines 
the error region where it crosses the unweighted line. 

which suggests that the CO gas is not at densities much 
greater than lO"*^* cm~^, or at temperatures much greater 
than 1000 K, though our models of the {v - 0) bands 
alone do not provide strong constraints on the high end 
of the temperature range. We find similar results within 
the errors for AA T au to the absorption line analysis of 
[France et all ()2012a '). Our 6- values are much higher than 
those found with sub-millimeter o bservations of the dis k 
of TW Hya (b < 0.04 km s"!) bv lHuehes et~all (|20T1 . 
however, this could be a consequence of their observa- 
tions probing colder, more quiescent CO gas at disk radii 
of ~ 100 AU. 

The errors are driven mainly by three sources: con- 
tinuum determination, low S/N, and saturation in the 
absorption lines. Estimates of the continuum fluxes are 
difficult due to the large number of emission lines in these 
spectra, especially the emission coincident with the ab- 
sorption bands. The error in the determination of the 
continuum flux is ^ 20%, leading to an error on the col- 
umn densities of ^ 0.25 dex. The low S/N in the contin- 
uum (S/N per resolution element ~ 7 for AA Tau) and 
absorption lines also contribute to the errors. Typical 
continuum flux values at A1420 A are listed in Table 3. 
At these low fluxes, the absorbers could not be detected 
with the Space Telescope Imaging Spectrograph (STIS), 
and have only recently been detected with COS. Satura- 
tion in these lines also leads to a degeneracy between the 
column densities and the Doppler &- values. We include 
higher-w bands in our fit as they have smaller oscillator 
strengths and are therefore less saturated. 

The line optical depths for all the targets (except HN 

Tau) are of order 100 - 600 for the \ow-v transitions {v 
= - 6) up to J-values of 15 to 25, which puts these lines 



on the flat portion of the COG. HN Tau, however, has 
optical depths of about 20 or less for all v transitions 
and even reaches optical depths of order unity in the 
models for the (7 - 0) band and therefore was undetected 
in the observations. The lower optical depths are caused 
by the higher temperature and higher fe- value of HN Tau. 
The higher temperature (700 K) of HN Tau lowers the 
line strength of each J-line compared to if the gas was 
cooler by distributing the absorption over more transi- 
tions. Similarly, the higher &-value (1.4 km s~^) low- 
ers the line center cross-section through larger Doppler 
broadening spreading out the cross-section in wavelength 
space. The product of the lower line strength and lower 
cross-section in each J-line leads to a significantly lower 
optical depth for HN Tau. None of the other targets ap- 
proach optical depths as small as 1 until the high- J levels 
of the V = 7, 8 states. 

4. RESULTS AND DISCUSSION 
4.1. CO Velocity and Isotopic Fraction 

The radial velocities of the CO absorption lines listed 
in Table 3 were obtained with a simple least squares fit. 
Assuming the best-fit parameters of the CO model, ve- 
locity shifts from -200 km s~^ to -1-200 km s"""^ with 1 
km s~^ intervals were explored. Due to low S/N, the (8 
- 0) band was not used in the velocity calculation. The 
average and standard deviation of the velocity shifts for 
the observed bands were taken as each target's velocity 
and velocity error. Shifting the velocity moves the model 
high- J CO lines relative to the observed spectra, which 
can change the best-fit temperature. However, the tem- 
perature would change by only ~ 50 K, which is well 
within our errors. The range of acceptable velocities was 
small enough that the best-fit parameters from the model 
were not affected by the shift. 

Since R W Aur contains str ong, redshifted H2 emis- 
sion lines (jFrance et al.ll2012bD . which sometimes places 
the emission on top of the CO absorption bandhead, we 
checked the velocity shift by hand. We required that the 
bandhead of the model and data match for the (1 - 0), 
(2 - 0), and (4 - 0) bands, which are the cleanest and 
strongest absorption bands. The best-fit velocity shift 
is +20 km s~^, consistent with th e stellar radial veloc - 
ity (-H15 km s~i) in the literature (|Nguven et al.ll2012D . 
The fit starts to become noticeably worse at ±20 km s~^ 
from the best-fit value, which we take as our error. 

The fitted CO absorption line velocities are generally 
consistent with the radial velocities of the stars from the 
literature to within the 15 km s~^ absolute uncertainty in 
the COS wavelength scale. However, the CO absorption 
in DE Tau appears to be somewhat blueshifted relative 
to the stellar velocity (see Tabl e 3), but is stil l consi stent 
within the errors. As noted bv lNguven et all (|2012[ ). the 
low stellar radial velocity of HN Tau (4.6 km s~^) devi- 
ates from the average velocity of the Taurus- Aurig a star- 
forming region (- 15 km s-i; IHartmann et al.| [l986'). 
The quoted stellar radial velocity of HN Tau is incon- 
sistent with our CO velocity, but is only slightly out- 
side our errors. These small velocity differences between 
the CO absorption and the star (Av = 2.04 - 23.4 km 
s~^) indicate that the CO is approximately at rest in the 
stellar frame, or at least not in a fast-moving disk wind. 
This is in agreement with studies presented bv lBast et al.l 
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(|2011f) and iNaiita et all (|2003D who find CO emission 
generally consistent with stellar velocities, with average 
velocity difference Av ~ 3.5 km s~^. However, the pos- 
sibility of a slow-moving disk wind, such as described in 
iPontop pidan et al.' (2011'), cannot be fully ruled out as 
an explanation for the location of the absorbing CO gas 
in our targets. 

The best-fit ^^CO/^'^CO ratio in our sample disks 
ranges from w 6 - 158 (See Table 4). The errors in 
these ratios, however, are very large because of the or- 
der of magnitude errors in the column densities. Our 
isotopic ratios are all con sistent with the in terstellar 
value of i^cO/i^CO ~ 70 (ISheffer et alllMl and the 
yo ung stellar obiect e nvironment value of ~ 100 found 
bv lSmith et~an ((20091 ) . but the S/N and spectral resolu- 
tion are not optimal for making a precise determination 
of the isotopic ratio in the disk. 

4.2. CO Temperature and Density 

We now compare our CO absorption temperatures and 
column densities to other gas measurements of the inner 
disk in the literature. The rotational excitation temper- 
atures of our disks (300 - 700 K) agree well with the 
analyses of U V fluorescent CO emission lines (460 ± 250 
K) in Schi ndhelm et all (|2012l ) for disks around similar 
CTTSs. In contrast, the average rotationa l temperature 
of the fuU disk sample of lSalvk et all (poTll ) is - 1100 K, 
albeit with a large range, which is higher than our sam- 
ple average of ~ 400 K. The difference is only marginally 
significant, although suggestive, due to the large un- 
certainties on the upper limit of the UV CO temper- 
ature distribution (see §3.3). For a direct comparison, 
the model fit to the AA Tau near-IR CO emission by 
I Salyk et al.l (|2011f) finds a rotation temperature of 950 K 
and logio(A^(CO))= 18.6, which are both hig her than our 
parameter values for the same target. iHorne et al.l ()2012f ) 
find a CO column density for A A Tau of 1.2x 10^^ cm~^ 
also h igher than our value. We compare the iSalvk et"al] 
(|2011f) and IHorne et all ()2012f ) CO parameters for A A 
Tau to our fit by using their column density and temper- 
ature values in our model and plotting them with ours 
in Figure 6. We assume a ^^CO/^'^CO ratio of 70, our b- 
value, and our velocity shift value for the other models. 
The other models do not appear to be consistent with 
our models or data. Because emission line studies probe 
a much larger region of the disk than the single line-of- 
sig ht of absorption lin e studies, it is not surprising that 
the ISalvk et al.l (|2011l ) study appears to probe a different 
molecular gas population. The higher temperature gas 
may be located at a disk height larger than the line-of- 
sight sampled in our absorption line measurements. The 
UV spectra may also be preferentially sensitive to lower 
column density gas due to extinction effects. However, 
the apparen t incon sistency between our model and the 
IHorne et all ()2012D model for AA Tau is surprising be- 
cause they are both absorption line studies. It may be 
that the near-IR continuum emission and the UV con- 
tinuum emission are produced in different locations, so 
that our lines-of-sight to them are different. Alterna- 
tively, if the CO is located in the inner disk, then the 
difference may be related to the warp of the inner disk, 
which varies with the rotation of the star ([Bouvier et al.l 
l2007f ). The temperature fits to our data show that it is 
inconsistent with a temperature of 2500 K, the H2 emis- 
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Fig . 6. — We use the t emperature a nd co lumn density values 
of Sal yk et al.l JMnT) and IHorne et all | |20TW) for AA Tau in our 
model and assume ^^CO/^^CO ~ 70 to compare their CO gas 
parameters to ours. For the comparison, we use the (4 - 0) band of 
A A T au. The data is in black, ou r fit is in orange, the ISalvk et al.l 
l|2011h model is in blue, and the IHorne et al.l l|2012l ) model is in 
red. 

sion/absorption line temperature found by France et al. 
(2012a, b) in these disks, and the H2 fiuorescence tem - 
per ature for TW Hya, model ed bv iHerczeg et al.l (|2004[ ) 
and INomura fc Millail ([2005). Assuming that both the 
CO and H2 are in local thermodynamic equilibrium, the 
hot H2 emission appears to be spatially separate from 
the CO absorption. 

Assuming that our CO rotational temperatures are 
representative of gas kinetic temperatures Tgas > 300 K, 
we estimate the maximum radius where warm CO can 
be maintained at this temperature (from Figure C.2 of 
iWoitke et al.l2011| ) to be rco ~ 10 AU This is consistent 
with t he UV CO emission studied bv I Schindhelm et"al] 
(|2012D who find that the emission arises from radii > 2 
AU in similar disks. Far-UV CO emission is spectrally 
unresolved from RECX-15, implying that rco > 3 AU. 
We estimate the number density of the CO in these disks 
by dividing our column densities by this radius: 



rco 

The estimated CO number densities are in the range of 
nco ~ 70 - 4000 cm~^, and are shown in Table 4. With 
these densities, we computed the CO/H2 ratio required 
for the CO to be thermalized up to J = 25 (the highest 
J-state that could be reliably identified in our data). If 
we assume that collisions with H2 are the leading con- 
tributor to the CO level populations, and that the H2 
density is sufficiently high, then the kinetic temperature 
will equal the CO rotational temperature, Trot{CO). Sig- 
nificant photoexcitiation would decouple the two temper- 
atures. We computed the H2 critical density (the den- 
sity where the spontaneous emission rate equals the colli- 
sional de-excitation rate) assuming an H2 ortho/para ra- 
tio of 3 and using collision rates calculated b y summing 
over all coUisional routes (|Yang et al.l 120101 ) downward 
out of level J = 25. The H2 critical density for CO ex- 
citation in the range 300 - 750 K is w (5.3 - 4.2) x 10^ 
cm~^. The CO/H2 ratios needed to thermalize the ab- 
sorption lines are listed in Table 4. They are upper limits 
and range from 1.6 x 10~^ — 7.9 x 10~^. For the CO den- 
sities derived with the assumed Woitke maximum radius, 



10 



McJunkin et al. 



DF Tau Non-Detection 



TABLE 4 

Protoplanetary Disk Warm Gas Parameters 



Object 


"CO 
cm~3 


CO/Ha^' 






AA Tau 


— 4.0 


< 1.1 X IQ- 


4 


r,t; 1 +373.0 

• -L 22 Q 


DE Tau 


i.i±;5:| X 


< 2.3 X IQ- 


4 


1 t- 0+778.5 
J-0.»_j4 2 


HN Tau 


6.6^26/ X IQl 


< 1.6 X IQ- 


5 


f5 0+93.7 
"■■^-5.0 


RECX-15 


4.2^3^8^ X 103 


< 7.9 X 10- 


4 


12 6+"3-3 


RW Aur 


5-3^4^6^ X 10^ 


< 9.9 X 10- 


5 


ICO c:+635.8 
iOO.O_j26.9 


SU Aur 


1.7tf;g X 103 


< 3.1 X 10- 


4 


.+2432.5 
'y-*_76.2 



^ CO/H2 ratio upper limit. For CO/H2 ratio loss than this value, the 
absorption lines arc thcrmalizcd up to J — 25. 

typical interstellar translucent and dense cloud CO/H2 
ratios (10"^ - 10'^: ILacv eraLlfTOOl iBurgh et al.|[2007t) 
are sufficient to maintain a thermal distribution for the 
observed CO absorption lines in five of the targets. Al- 
though the derived CO densities for most of our targets 
suggest that the CO is thermalized, the sample is also 
consistent with sub-thermal excitation, within the un- 
certainties. Sub-thermally populated lines would lead to 
the CO temperatures being underestimated. If the ab- 
sorbing CO population is in thermal equilibrium, then 
we conclude that the observed CO/H2 ratios derived by 
France et al. (2011a, 2012a) are not representative of the 
local CO/H2 ratios in the warm molecular disk surface. 

We compare our derived temperatures and approxi- 
mated CO densities (which we change to a molecular hy- 
drogen density with an assumed CO/H2 ratio of lO"**) 
to the Woitkc et al. (2011) model of RECX-15 as a check 
on the vertical disk structure. Continuing our assump- 
tion that collisions with H2 are the leading contributor to 
the CO level populations, the averaged hydrogen density 
will be dominated by H2 and we can directly compare 
the plots in Figure C.2 of iWoitke et all (|201ll ) with our 
temperatures and densities assuming the CO radius to be 
Tco ~ 10 AU as above. For each target's temperature 
and density pair, the height of the gas determined from 
the model are consistent with each other. These heights 
range between z/r ^ 0.6 to ~ 0.7 for the six targets, 
putting the gas in the flared upper disk atmosphere. 

4.3. Comparison of CO with System Parameters 

Using the scale height model from §4.2, we calculate 
that only inclinations of > 79° will intercept the Ay = 1 
surface of the disk, which is larger than the inclinations 
in Table 1. At the inclinations of the targets, the sight- 
lines are not probing into the depth of the disk where 
the visual extinction exceeds unity, which indicates that 
our CO is likely well above the Ay = 1 disk surface. 
However, using a dimensionless height oi z/r — 0.6, the 
disk would be intercepted for inclinations > 60°. This 
increases our estimate on the unconstrained HN Tau in- 
clination. The lower limit on inclinations of > 60° is com- 
patible with most of our disk inclinations, with the no- 
table exception of DE Tau. The inclination of DE Tau is 
too small for us to be observing the gas in absorption, yet 
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Fig. 7.— COS data for DF Tau with the model fit of AA Tau 
overplotted, showing the absence of CO absorption in the DF Tau 
spectrum (DF Tau data in black and best-fit model for AA Tau 
m orange ). The "CO bandheads are marked to show the two 
different CO species. 



we clearly see the absorption bands in the data. DE Tau 
may be a candidate for a protoplanetary disk with a slow- 
moving molecular disk wi nd with a more face-on in clina- 
tion and blue-shifted lines (iPontoppidan et al.ll2Gllf). DF 
Tau, another target in the France et al. "(i2012bl s ample 
known to have CO (iNai ita et al.ll2003l: ISalvk et al.ii2011l: 



ISchindhelm et al.ll2012l) and thought to h ave a high inch 
nation f85°: iJohns-KruU fc Valentill2001[ ). does not show 
UV CO absorption, which is surprising. DF Tau has an 
anomalously large ra dius for a 10^ Myr star (3.37 i?©, 
iGuUbring et al.lfl998| ). which leads to the large inclina- 
tion value. The spectrum for DF Tau is shown in Fig- 
ure 7 with the fit for AA Tau overplotted to illustrate 
the non-detection. A revised determination of the incli- 
nation of DE Tau and DF Tau, and possibly the other 
sources in the sample, would be useful. With notable ex- 
ceptions, for our mid- to high-inclination disks known to 
contain CO, we see absorption if the gas is located at a 
height of z/r > 0.6, which is compatible with our derived 
temperatures and estimated densities. A cartoon of this 
geometry is shown in Figure 8. 

In Figure 9, we compare the column density and rota- 
tional temperature of CO with the mass accretion rate 
and disk inclination. We do not see any strong correla- 
tions between these parameters, though there may be a 
slight decrease in temperature with increasing disk incli- 
nation. At higher inclinations, we are observing the disk 
more edge-on so that the line-of-sight passes through the 
colder, denser material lower in the disk atmosphere. 

5. CONCLUSIONS 

We present model fits to ultraviolet CO absorption 
lines in six protoplanetary systems. We find results in 
AA Tau and HN Tau that are consistent with previ- 
ous work and extend our analysis to four new targets 
obtained with HST-COS far-UV G130M and G160M 
modes. We find CO rotational temperatures in the range 
300 - 700 K, which agree well with UV CO fluorescence 
rotational temperatures. However, our temperatures are 
cooler than inner disk CO temperatures obtained from 
the analyses of IR CO e mission spectra and warmer than 
the iRettig et all ()2006D CO temperatures from IR CO 
absorption. The IR emission data may be probing a dif- 
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Fig. 8. — A schematic representation of the line-of-sight geome- 
try for the inner region of T Tauri star disks. The red layer is ob- 
served as low-density , hot (T ~ 2500 K) H2 emission and absor ption 
IIHerczeg et al.l [2001 : lYang et all [20Tll : [France et all I2012al '): the 
green material is observed as intermediate density, warm (T ~ 500 
K) CO emission and absorption; and the dark midplane is inacces- 
sible to UV observations. 
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ferent gas population at smaller disk radii, while the IR 
absorption data is likely probing gas at larger radii. The 
measured velocities of the CO absorbing gas rule out an 
origin in a fast-moving disk wind. Our derived temper- 
atures and approximated densities of the gas are con- 
sistent with models for disk heights of zjr ^ 0.6. This 
CO location constrains the inclination of our disks to 
be > 60° in order to intercept the absorbing gas in the 
model. 

We note that the present analysis is roughly at the 
limit of what is feasible with the current generation of 
HST instrumentation. Higher spectral resolution would 
improve our molecular parameter determination signifi- 
cantly. Unfortunately, observations at these flux levels 
with HST-STIS E140M mode are not feasible. A new 
observational capability will be required to derive more 
robust disk parameters from CO absorption line obser- 
vations of CTTSs. However, future observations with 
the new COS G130M A1222 mode should be able to di- 
rectly probe warm H2 absorption in these disks at wave- 
lengths around 1100 A. These observations would allow 
us to directly determine the value of the CO/H2 ratio in 
moderate-to-high inclination protoplanetary disks. 

We thank M. Eidelsberg for providing the CO oscilla- 
tor strengths and rovibrational line wavelengths used in 
this work. This work made use of data from HST guest 
observing program 11616 and was supported by NASA 
grants NNX08AC146 and NAS5-98043 to the University 
of Colorado at Boulder. 
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Correlations of the target parameters and best-fit pa- 



REFERENCES 



Akeson, R. L., Ciardi, D. R., van Belle, G. T., & Creech-Eakman, 

M. J. 2002, ApJ, 566, 1124 
Andrews, S. M. &; Williams, J. P. 2007, ApJ, 659, 705 
Armitage, P. J., Clarke, C. J., & Palla, F. 2003, MNRAS, 342, 

1139 

Bary, J. S., Weintraub, D. A., & Kastner, J. H. 2003, ApJ, 586, 
1136 

Bary, J. S., Weintraub, D. A., Shukla, S. J., Leisenring, J. M., &; 
Kastner, J. H. 2008, ApJ, 678, 1088 

Bast, J. E., Brown, J. M., Herczeg, G. J., van Dishoeck, E. F., & 
Pontoppidan, K. M. 2011, A&A, 527, A119 

Bergin, E., Calvet, N., Sitko, M. L., Abgrall, H., D'Alessio, P., 
Herczeg, G. J., Roueff, E., Qi, C, Lynch, D. K., Russell, R. W., 
Brafford, S. M., & Perry, R. B. 2004, ApJ, 614, L133 

Bertout, C, Basri, G., & Bouvier, J. 1988, ApJ, 330, 350 

Bertout, C, Robichon, N., & Arenou, F. 1999, A&A, 352, 574 

Bitner, M. A., Richter, M. J., Lacy, J. H., Herczeg, G. J., 
Greathouse, T. K., Jaffe, D. T., Salyk, C, Blake, G. A., 
HoUenbach, D. J., Doppmann, G. W., Najita, J. R., & Currie, 
T. 2008, ApJ, 688, 1326 

Bouvier, J., Alencar, S. H. P., Boutelier, T., Dougados, C, Balog, 
Z., Grankin, K., Hodgkin, S. T., Ibrahimov, M. A., Kun, M., 
Magakian, T. Y., & Pinte, C. 2007, A&A, 463, 1017 



Bouvier, J., Grankin, K. N., Alencar, S. H. P., Dougados, C., 
Fernandez, M., Basri, G., Batalha, C, Guenther, E., Ibrahimov, 
M. A., Magakian, T. Y., Melnikov, S. Y., Petrov, P. P., Rud, 
M. v., & Zapatero Osorio, M. R. 2003, A&A, 409, 169 

Burgh, E. B., France, K., & McCandliss, S. R. 2007, ApJ, 658, 446 

Calvet, N. & GuUbring, E. 1998, ApJ, 509, 802 

Carmona, A., van den Ancker, M. E., Henning, T., 

Pavlyuchenkov, Y., DuUemond, C. P., Goto, M., Thi, W. F., 
Bouwman, J., & Waters, L. B. F. M. 2008, A&A, 477, 839 

Carpenter, K. G., Robinson, R. D., Wahlgren, G. M., Linsky, 
J. L., & Brown, A. 1994, ApJ, 428, 329 

Carr, J. S., Tokunaga, A. T., Najita, J., Shu, F. H., & Glassgold, 
A. E. 1993, ApJ, 411, L37 

Danforth, C. W., Keeney, B. A., Stocke, J. T., ShuU, J. M., & 
Yao, Y. 2010, ApJ, 720, 976 

Donati, J.-F., Skelly, M. B., Bouvier, J., Gregory, S. G., Grankin, 
K. N., Jardine, M. M., Hussain, G. A. J., Menard, F., 
Dougados, C, Unruh, Y., Mohanty, S., Auriere, M., Morin, J., 
Fares, R., & MAPP Collaboration. 2010, MNRAS, 409, 1347 

Durrance, S. T. 1981, J. Geophys. Res., 86, 9115 

Eidelsberg, M., Jolly, A., Lemaire, J. L., Tchang-Brillet, W.-t)., 
Breton, J., & Rostas, F. 1999, A&A, 346, 705 

Eidelsberg, M. & Rostas, F. 2003, ApJS, 145, 89 

Eisner, J. A., Hillenbrand, L. A., White, R. J., Bloom, J. S., 
Akeson, R. L., & Blake, C. H. 2007, ApJ, 669, 1072 



12 



McJunkin et al. 



Fedele, D., van den Anckcr, M. E., Hcnning, T., Jayawardhana, 

R., & Olivcira, J. M. 2010, A&A, 510, A72 
Federman, S. R., Glassgold, A. E., Jenkins, E. B., & Shaya, E. J. 

1980, ApJ, 242, 545 
Feldman, P. D. & Brune, W. H. 1976, ApJ, 209, L45 
Feldman, P. D., Burgh, E. B., Durrance, S. T., & Davidsen, A. F. 

2000, ApJ, 538, 395 
France, K., Burgh, E. B., Herczeg, G. J., Schindhelm, E., Yang, 

H., Abgrall, H., Roueff, E., Brown, A., Brown, J. M., & Linsky, 

J. L. 2012a, ApJ, 744, 22 
France, K., Schindhelm, E., Burgh, E. B., Herczeg, G. J., Harper, 

G. M., Brown, A., Green, J. C., Linsky, J. L., Yang, H., 
Abgrall, H., Ardila, D. R., Bergin, E., Bethell, T., Brown, 

J. M., Calvet, N., Espaillat, C, Gregory, S. G., Hillenbrand, 
L. A., Hussain, G., Ingleby, L., Johns-KruII, C. M., RouefT, E., 
Valenti, J. A., & Walter, F. M. 2011a, ApJ, 734, 31 
France, K., Schindhelm, E., Herczeg, G. J., Brown, A., Abgrall, 

H. , Alexander, R. D., Bergin, E. A., Brown, J. M., Linsky, 
J. L., Roueff, E., & Yang, H. 2012b, ApJ, 756, 171 

France, K., Yang, H., & Linsky, J. L. 2011b, ApJ, 729, 7 

Froning, C. S., Cantrell, A. G., Maccarone, T. J., France, K., 
Khargharia, J., Winter, L. M., Robinson, E. L., Hynes, R. I., 
Broderick, J. W.. Markoff, S., Torres, M. A. P., Garcia, M., 
Bailyn, C. D., Prochaska, J. X., Werk, J., Tliom, C., Bcland, S., 
Danforth, C. W., Keeney, B., & Green, J. C. 2011, ApJ, 743, 26 

Green, J. C., Froning, C. S., Osterman, S., Ebbcts, D., Heap, 
S. H., Leitherer, C., Linsky, J. L., Savage, B. D., Sembach, K., 
ShuUj J. M., Siegmund, O. H. W., Snow, T. P., Spencer, J., 
Stern, S. A., Stocke, J., Welsh, B., Bcland, S., Burgh, E. B., 
Danforth, C., France, K., Keeney, B., McPhate, J., Penton, 
S. v., Andrews, J., Brownsberger, K., Morse, J., & Wilkinson, 
E. 2012, ApJ, 744, 60 

Gullbring, E., Calvet, N., MuzeroUe, J., & Hartmann, L. 2000, 
ApJ, 544, 927 

Gullbring, E., Hartmann, L., Briceno, C., & Calvet, N. 1998, 

ApJ, 492, 323 
Haridass, C. & Huber, K. P. 1994, ApJ, 420, 433 
Hartigan, P., Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736 
Hartmann, L., Hewett, R., Stabler, S., & Mathieu, R. D. 1986, 

ApJ, 309, 275 

Herczeg, G. J., Wood, B. E., Linsky, J. L., Valenti, J. A., & 

Johns-KruU, C. M. 2004, ApJ, 607, 369 
Hernandez, J., Calvet, N., Briceno, C, Hartmann, L., Vivas, 

A. K., MuzeroUe, J., Downes, J., Allen, L., & Gutermuth, R. 

2007, ApJ, 671, 1784 
Home, D., Gibb, E., Rettig, T. W., Brittain, S., Tilley, D., & 

Balsara, D. 2012, ApJ, 754, 64 
Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2005, Icarus, 

179, 415 

Hughes, A. M., Wilner, D. J., Andrews, S. M., Qi, C., & 

Hogerheijde, M. R. 2011, ApJ, 727, 85 
Ida, S. & Lin, D. N. C. 2004, ApJ, 604, 388 
Ingleby, L., Calvet, N., Bergin, E., Herczeg, G., Brown, A., 

Alexander, R., Edwards, S., Espaillat, C, France, K., Gregory, 
S. G., Hillenbrand, L., Roueff, E., Valenti, J., Walter, F., 
Johns-KruU, C, Brown, J., Linsky, J., McClure, M., Ardila, D., 
Abgrall, H., Bethell, T., Hussain, G., & Yang, H. 2011, ApJ, 
743, 105 

Johns-KruU, C. M. & Valenti, J. A. 2001, ApJ, 561, 1060 
Johns-KruU, C. M., Valenti, J. A., & Linsky, J. L. 2000, ApJ, 
539, 815 



Kraus, A. L. & Hillenbrand, L. A. 2009, ApJ, 704, 531 

Lacy, J. H., Knacke, R., Geballe, T. R., & Tokunaga, A. T. 1994, 

ApJ, 428, L69 

Lawson, W. A., Lyo, A.-R., & MuzeroUe, J. 2004, MNRAS, 351, 
L39 

Loinard, L., Torres, R. M., Mioduszewski, A. J., Rodriguez, L. F., 
Gonzalez-Lopezlira, R. A., Lachaume, R., Vazquez, V., & 
Gonzalez, E. 2007, ApJ, 671, 546 
Lopez-Martin, L., Cabrit, S., & Dougados, C. 2003, A&A, 405, LI 
Luhman, K. L. & Steeghs, D. 2004, ApJ, 609, 917 
Mamajek, E. E., Lawson, W. A., & Feigelson, E. D. 1999, PASA, 
16, 257 

McPhate, J. B., McCandliss, S. R., Feldman, P. D., Burgh, E. B., 
& Pelton, R. 1997, in Bulletin of the American Astronomical 
Society, Vol. 29, AAS/Division for Planetary Sciences Meeting 
Abstracts #29, 1051 
Najita, J., Carr, J. S., & Mathieu, R. D. 2003, ApJ, 589, 931 
Nguyen, D. C, Brandeker, A., van Kerkwijk, M. H., & 

Jayawardhana, R. 2012, ApJ, 745, 119 
Nomura, H. & Millar, T. J. 2005, A&A, 438, 923 
Pascucci, I., Gorti, U., HoUenbach, D., Najita, J., Meyer, M. R., 
Carpenter, J. M., Hillenbrand, L. A., Herczeg, G. J., Padgett, 
D. L., Mamajek, E. E., Silverstone, M. D., Schlingman, W. M., 
Kim, J. S., Stobic, E. B., Bouwman, J., Wolf, S., Rodmann, J., 
Hincs. D. C, Lunine, J., & Malhotra, R. 2006, ApJ, 651, 1177 
Pontoppidan, K. M., Blake, G. A., & Smette, A. 2011, ApJ, 733, 
84 

Pontoppidan, K. M., Blake, G. A., van Dishoeck, E. F., Smette, 
A., Ireland, M. J., & Brown, J. 2008, ApJ, 684, 1323 

Ramsay Howat, S. K. & Greaves, J. S. 2007, MNRAS, 379, 1658 

Rettig, T., Brittain, S., Simon, T., Gibb, E., Balsara, D. S., 
Tilley, D. A., & Kulesa, C. 2006, ApJ, 646, 342 

Ricci, L., Testi, L., Natta, A., Neri, R., Cabrit, S., & Herczeg, 

G. J. 2010, A&A, 512, A15 

Salyk, C, Blake, G. A., Boogert, A. C. A., & Brown, J. M. 2011, 
ApJ, 743, 112 

Schindhelm, E., France, K., Burgh, E. B., Herczeg, G. J., Green, 
J. C, Brown, A., Brown, J. M., & Valenti, J. A. 2012. ApJ, 
746, 97 

Sheffer, Y., Rogers, M., Federman, S. R., Lambert, D. L., & 

Gredel, R. 2007, ApJ, 667, 1002 
Smith, R. L., Pontoppidan, K. M., Young, E. D., Morris, M. R., 

& van Dishoeck, E. F. 2009, ApJ, 701, 163 
Trilling, D. E., Lunine, J. I., & Benz, W. 2002, A&A, 394, 241 
Vidal-Madjar, A., Lagrange-Henri, A.-M., Feldman, P. D., Beust, 

H. , Lissauer, J. J., Deleuil, M., Ferlet, R., Cry, C, Hobbs, 

L. M., McGrath, M. A., McPhate, J. B., & Moos, H. W. 1994, 
A&A, 290, 245 

White, R. J. & Ghez, A. M. 2001, ApJ, 556, 265 

Woitke, P., Riaa, B., Duchene, G., Pascucci, I., Lyo, A.-R., Dent, 
W. R. F., Phillips, N., Thi, W.-F., Menard, F., Herczeg, G. J., 
Bergin, E., Brown, A., Mora, A., Kamp, I., Aresu, G., Brittain, 
S., de Gregorio-Monsalvo, I., & SandeU, G. 2011, A&A, 534, 
A44 

Yang, B., Stancil, P. C, Balakrishnan, N., & Forrey, R. C. 2010, 

ApJ, 718, 1062 
Yang, H., Linsky, J. L., & France, K. 2011, ApJ, 730, LIO 



